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The eta-photon transition form factor is evaluated in a formalism based on a phenomenological 
description at low values of the photon virtuality, and a QCD-based description at high photon 
virtualities, matching at a scale Qq. The high photon virtuality description makes use of a Distribu- 
tion Amplitude calculated in the Nambu-Jona-Lasinio model with Pauli-Villars regularization at the 
matching scale Qq, and QCD evolution from Qq to higher values of Q 2 . A good description of the 
£Sj . available data is obtained. The analysis indicates that the recent data from the BaBar collaboration 

on pion and eta transition form factor can be well reproduced, if a small contribution of higher twist 
is added to the dominant twist two contribution at the matching scale Qq. 
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I. INTRODUCTION 



J2' 

Meson Distribution Amplitudes (DA) are fundamental theoretical ingredients in the description of exclusive high 
JL . energy processes. The pseudoscalar transition form factors (TFF), F 77 *p, describing the process P — > 77*, where P 
is a pseudoscalar meson, are directly connected with the DAs. Recently, the BaBar Collaboration has provided new 
data at high virtuality for the pion and eta TFF (irTFF and r]TFF) 0, Q ■ The implications of these results on 
'— irTFF in our understanding of the pion structure have been widely discussed [H-Qjl 

In particular, these results have cast doubts on the behavior, as a function of the light-cone momentum fraction x, 
of the pion distribution amplitude (ttDA) 4br(x) [3 EH, a quantity for which some investigations have predicted a flat 
^vq I behavior, i.e., a constant value for any x [J, m good agreement with the data of the form factor. These scenarios 
are compatible with QCD sum rules [l5[ and lattice QCD 0, [ItJ calculations which provide values for the second 
moment of the itDA which are large compared to the asymptotic value 6x(l — x). Several model calculations, such 
^*>0 as the ones performed in the Nambu-Jona-Lasinio (NJL) (l8rl2fj| or in the "spectral" quark model [2l[ frameworks, 

■ give a constant nDA, i.e. ^(x) = 1. 
t-H \ With the availability of data about the eta, it is important to analyze all the proposed theoretical schemes. Some 
work in this direction has been already done [TJ l23 - l27j . In particular, in Refs. [ill, I23 - |25| the importance of the 
transverse momentum of the quarks is emphasized, making use of some parametrization of the eta wave function. 
In Refs. [2(1 [27|, the TFFs of pseudoscalar mesons are considered by using a dispersive representation of the axial 
anomaly, considering also the violation of factorization and possible higher twist corrections. 

The parton distributions, generalized parton distributions and distribution amplitudes have been used as a test of 
hadron models. The procedure consists of three ingredients: i) the hadron model provides a low energy description 
of the studied distribution; ii) a high energy description is obtained by QCD evolution, which needs an input at 
some low scale Qq; Hi) a matching condition between the two descriptions at a scale Qq characterizing the separation 
between the two regimes. This pro cedure has been useful in the study of nucleon parton distributions |28|-|30| as well 
as in that of pion distributions |3ll43a |. 

In Ref. 0, a version of the previous program, but in a rather model- independent formalism, has been used to 
calculate the itTFF. An excellent description of experimental data has been obtained in the whole range of virtuality. 
Summarizing, the evaluation of the irTFF at high Q 2 values in Ref. [8j is based on the following arguments: 
i) chiral symmetry and soft pion theorems, which explain that, at some point Qq, the irDA has a flat behavior, 
4>Tr(x, Qq) = 1; ii) applying QCD evolution to the ttDA, one can obtain the irTFF at any Q 2 ^> Qq, Hi) for Q 2 < Qq, 
the experimental parametrization of i 7 77 * w (Q 2 ) given in ref. [36| is assumed; iv) for Q 2 > Qq the ttTFF is given 
by its standard expression in terms of the ttDA modified in two directions, the quark propagator is corrected , as 
suggested by Radyushkin Q, and a term originated by other higher twist contributions is included. 
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This scheme, successful for the pion, requires further tests. The most natural one consists in the evaluation of 
the same quantity for the other pseudoscalar mesons. In this paper, the program is developed for the 77 meson, for 
which a few sets of data are available [H, [3(| ■ The importance of the 77 — 77' system for our understanding^ the 
QCD symmetries, and for their treatment in effective, low energy descriptions, is well known (see, i.e., Ref. [37[ and 
references therein). To implement this program, the approach of Ref. [8[ has to be complicated, and some hints 
have to be obtained within a specific model. In particular, a generalized SU(3) Nambu-Jona-Lasinio model with 
Pauli-Villars regularization, along the lines of Ref. [38[ , will be used. 

The NJL model is the most realistic model for the pseudoscalar mesons based on a local quantum field theory built 
with quarks. It respects the realization of chiral symmetry and gives a good description of the low energy physics of 
pseudoscalar mesons. It allows to describe mesons in a field theoretical framework treating them as bound states in 
a fully covariant manner using the Bethe-Salpeter amplitude. In this way, the Lorentz covariance of the problem is 
preserved. 

The NJL model is a non-renormalizable field theory and therefore a cut-off procedure has to be implemented. The 
Pauli-Villars regularization procedure has been chosen because it respects the gauge symmetry of the problem. The 
NJL model together with its regularization procedure is regarded as an effective theory of QCD. In the chiral limit, 
it predicts for the pion <f>^ (x) = l,0<x<l,in agreement with the model independent study of the pion in Ref. 0. 
At this point, due to the lack of fits in the SU(3) NJL model with the Pauli-Villars regularization , it is performed a 
new analysis of the parameters of the model. It is interesting to notice that an early attempt to use the NJL model, 
but in a L^(3) invariant scheme, has been presented in Ref. [39j, where the 77 parton distribution has been evaluated. 

The paper is organized as follows. In section II, the theoretical description of the 77* — > P is reviewed, extracting 
the soft (non-perturbative) part, to be described by the model. In section III, the 77 calculation in the NJL model is 
presented. In the following section, numerical results are presented and discussed. The conclusions are drawn in the 
last section. In the Appendix, a summary of the NJL model is given, including the description of the fit which has 
been used. 



II. THE 77* -> P PROCESS: THEORETICAL DESCRIPTION. 



The subject of this study is the transition form factor, _F 77 »p, i.e., the form factor for the coupling of a real photon 
and a virtual photon to a pseudoscalar meson, P. The TFF is a very important quantity in the QCD description of 
exclusive processes. In particular, it can be used to obtain information on the shape of the meson DA [Til liol. l4l|. 
Experimentally, it has been measured for the 7r, 77 and 77' mesons by the CELLO by the CLEO [36[ and, recently, 
by the BaBar [l|, collaborations. The latter results, for the pion, have been found in disagreement with theoretical 
expectations. 

In order to establish the proper formalism, in this section the theoretical description of the 7 (q\) 7* ((72) P {k) 
process is reviewed. From general arguments it is well known that the transition amplitude of this process can be 
written as (see, i.e., [Hj]): 

(P (k) out| 7 { gi , e x ) , 7* ((72, £2) in) = 1 (2tt) 4 S 4 {q x + q 2 - k) T , (1) 

with 



T = iTrae lf ,e 2l/ q lp q 2 ae flupa F^p , (2) 

where a is the fine structure constant. On the other hand, applying the reduction formalism of Lehmann, Symanzik 
and Zimmermann 43], this process is described by 

(P(fc)out|7(gi,ei),7* (q 2 ,£ 2 )m) 

= -(2nf5 i (q 1 +q 2 -k) J d 4 z e "*(«-*)*/a {P (fc) out | T ( £l . j (£) £2 . j (_£) ) | ) , (3) 

where (z) = ej ipj (z) (z) is the electromagnetic current for the quarks. Evaluating the time ordered 
product at leading order and after some algebra, one has 



dH 



e -itz 



]T 2 i e] ElM e 2v (-ie^H a ) f2 _ ^ + ie < P ( fc ) out l ^ (f ) W^j (™) |0> • (4) 
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P(fc)out|^ a (|) (-|) |0) = f -^e-^-^i^il) , (5) 



Defining 

dH 

where a, b are quadrispinor indices, one can write the amplitude T as, 

T = e^e 2lJ e^ I a p , (6) 

with 

I aP = -i f 2 («!-/)„ J e * (^\ a (^i b (i)) ■ (7) 

J (2tt) (qi -t) -mq + ie V ' 

Since the symmetric part of / a/ g doesn't give contribution to 7", the attention is focused in the antisymmetric part. 
With the momenta q\ and q 2 , one can build two antisymmetric tensors, q\ a q 2 p — qip q 2a and £ a ppa-qiq 2 ■ Nevertheless, 
there is not enough structure in the integrand of Eq. ([7]) to generate a tensor like £ a ppa<li q 2 , at least at the leading 
order. Therefore, the tensor structure of I a p is 

£*/3 = ijj (qia 92/3 - qip Q2a) I + symmetric part . (8) 

Turning back to the Eqs. © and ([5]), and using Eq. ©, one gets F 77 »p = —I. Contracting flSJ) with q 2 p — qxp q 2a 
and using the explicit expression of I a p, Eq. (O, yields 



f' ; 



77 * P (qx.^) 2 7 (2TT) 4 Y - - + l€ ^ 

lqi-(qi-e)(42j 5 ) ba -q2.(qi-e)(4ij 5 ) ba \(txi b ( £ )) • ( 9 ) 

This expression for the transition form factor is quite general. At this stage the assumptions made are i) the free 
quark propagator has been used in going from Eq. (J3|) to Eq. (U) and ii) the corrections to the elctromagnetic vertex 
has not been considered. A more general expression could be obtained by changing the free propagator, (f — m + ie) 1 , 
by the general one associated to a dressed quark, (A (t) f — B (t) + ie) , studied in actual lattice QCD calculations 
44] and including a term with the neglected structure, s a p p(y qiq 2 . 

Looking at the kinematics of the process, one can choose the reference frame in such a way that the pion and 
photons four momenta arc k = (Ek, k x , 0, k z ) , q\ — (E%, 0, 0, — E\) , q 2 = (Ek — E\, k x , 0, E\ + k z ) , respectively. 
It is interesting to express the quantities in terms of the light-front variables, fc ± = (Ek ± k z ) /V%, k± — (k x ,0) , 
<h = {Q 2 + m p) / 2fc+ > qt = °> 5i-L =0, q£ = k + , q 2 =- (Q 2 - k\) /2fc+, q 2± = k±, where Q 2 = -q\. In the limit 
of large Q 2 , some of the quantities in Eq. ([9]) can be approximated by 



ti 



Q 



2 



7 + , (10) 



qi .q 2 ~ X -Q 2 , (11) 

( qi -£) 2 - m 2 + tec-^-Q 2 , (12) 



giving for the transition form factor the expression 

1 f d*£ ^ 1 o / 1 



Finally, defining x = one arrives to the usual expression 



) (V nb W) ■ (13) 

/ ha 



1 /"da; 

Fy^*p(Q ) ~ ^2 / — <J>p(x,Q 2 ) , (14) 
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with 



= -2i — — e 



(27T) 

«-* + (-*) (P(*)out|$>2& (I) 7 + 75 ^- (-|) |0) 



2tt 



(15) 



As will be discussed later, in Eq. (fT4"|) . besides the explicit Q 2 dependence, also an implicit one appears, through the 
QCD evolution ol $p (a;, Q 2 )- In the SU(3) formalism, the quark operator has the form 

E^i (|) (-f) = ^ (|) 7 + 75^ (-§) , (16) 

i 

with = 2 A°/(3v / 6) + A 3 /6 + A 8 /(6v / 3) where A a are the SU(3) generators. In the present case it is more interesting to 
use the flavor basis in describing the rj particle (see the Appendix). In this basis, one has O = 5 A g /18+A s / (9v / 2)+A 3 /6. 
As usual, the DA of P in the flavor basis is defined as 



l f P ^Q 2 ) = -J ^e**- fc+ (*-*) (P(fe)|^-(f)7 + 7b^(-f)|0> 
with j = 3, q, s. This yields 



(17) 

z+=0,z±=0 



<S> P (x, Q 2 ) = ^fU%& Q 2 ) + |V2/^^(x, Q 2 ) + -fMUx, Q 2 ) . (18) 



In the pion case, this equation corresponds to <£v (x) = y2 /„• (/>„• (a;) /3, where </> ff (a;) is the 7T-D^4 and /„■= 0.131 GeV 
is the pion decay constant. 

One should notice that, in going from Eq. ([!} to the final result Eq. (]14p. a few approximations have been done: the 
free expression has been used for the quark propagator, with the additional simplification given by Eq. (|12p ; besides, 
the approximations Eqs. (|10[) and have been applied in the numerator of Eq. ([5]) and a new tensor structure in 
I a p has been neglected. Some of these corrections have a kinematic character, while others are certainly dynamical. 
Both type of corrections imply the presence of higher twist distribution amplitudes. 

In Ref. [8| it has been argued that the approximations, leading to the simple expression Eq. (|14p for the transition 
form factor, are too crude to explain the BaBar experimental data, and corrections at the next order in the Q~ 2 
expansion have been added. The simplest way to implement these corrections is to start from the following expression: 



Q 2 F 77 , P (Q 2 ) = / ^t$p(x, Q 2 ) + % . (19) 

X + V 



The mass M in Eq. ([15)1 was introduced by Radyushkin Q, to cure the divergence of the integrand in Eq. (fT4l . 
occurring when a DA $p (x,Qq), not vanishing at x = 0,1, is used. This was justified as a consequence of the 
existence of some transverse component in the quark momentum. As it has been shown in the previous section, M 
contains not only effects associated to the mean transverse momentum, but also the ones associated to the constituent 
quark masses, among others. In Ref. [§[ it has been shown that it is necessary to introduce the C3-dependent term in 
Eq. (|19p. otherwise the data cannot be well described around the region of Q 2 = 10- 20 GeV 2 . The inclusion of this 
term has been thoroughly motivated in this section, where it has been shown that the perturbative approach leading 
to Eq. (|14p is correct only for high enough values of the virtuality. We call the C3 term as "the higher twist term" , 
although it is clear that also the mass term, M, is of the same order. 

III. THE r) TRANSITION FORM FACTOR 

In this section, we evaluate the TjTFF. To this aim, according to Eq. ([19]). in calculating F 11 * rj (Q 2 ) the rjDA is 
needed. From Eq.(|18p. the r]D A is expressed by 

% (x, Q 2 ) = ^V2f* 0« (x, Q 2 ) + 2 - f v s (x, Q 2 ) . (20) 
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In the calculation, we use the following values of the rj weak decay constants 

/« = (0.828 ±0.019)/, , 

/' = (-0.848 ±0.042)/,, (21) 

with f n = 131 MeV, obtained in the phenomenological study of Ref. [37j . 

Now, it is necessary to calculate the tjDA at some initial scale Qo within a model. To this end, we have obtained 
the DAs corresponding to the q and s flavors within the Nambu-Jona Lasinio (NJL) model, which has a long tradition 
of successful predictions of meson parton structure (3l| - [35j . In particular, we use in the present calculation the three 
quark flavor version of the model with the Pauli-Villars regularization [38l |45( . A brief summary of the model and of 
the regularization procedure is given in the Appendix. 

In the NJL model, mesons are described through Bethe-Salpeter amplitudes. For the rj meson one has: 

(77 (k)\ff>p (a*) (a*) |0) = e **C*i+*»V 2 J ALj^-**) 



iS r ( I - ^ j ig mq {cos4>\ q - sin 4> X s ) ij 5 iS F ( £ + 7; 



(22) 

Pa 



where the index a, j3 stand for color, flavour and quadrispinor index. Inserting this expression in Eq. (JTTJ) one obtains 
if* (x) = V2g vqq (cos - smc/)8j tS ) N c 

<fM' + ~ * + {' -l)) tT {'-I) ■>•* (' + \ ) ^) ■ <23) 

Evaluating the trace and using Eq. (|A18[) of the Appendix, one has 

<t> q n {s) (x, Qo) = —, 2T ( x > m <?W ' m %) ' ( 24 ) 

i 2 (m g ( s ),m g ( s ),m^) 

where I2 (x,mi,rnp) is given in Eq. (|A23I) and l2( m q(s) j ^ l q(s) i ^«) * s the two-propagator integral defined in Eq. 
V2[) . The flavor I? A defined in Eq. ([24]) satisfies the normalization condition: 

^f)( S ,Q2) = i. (25) 



For Q 2 > Q 2 ,, the rjTFF is obtained through QCD evolution [Tj, El, El- Tne can be expressed in terms of 

the Gegenbauer polynomials, 



n(even)—0 



log 



Q 2 



' s (x,Q 2 ) = 6x(l-x) Yl < s C 3 J 2 (2x-l)\ , (26) 



where 7„ is the anomalous dimension 



°F (, A ^ 1 



(n + l)(n + 2) / ' 



(27) 



/3 = 11 ^ c ' — is the beta function to lowest order and Cf — If the a^f coefficients are known, using Eq. 

(|26|) in Eq. (|20|) and (fl9|) . the rjTFF is obtained for any Q 2 3> Aq Cd . Once 4>% s (x, Q 2 ,) are known at a given scale 
Qo, the a*' s coefficients are obtained using the orthogonality property of the Gegenbauer polynomials 



2n + 3 



3 (n + 1) (n + 2) J 



dxC 3 J 2 (2x-l)^ s (x,Q 2 ) . (28) 



In this scheme, the 77 meson cannot couple to two gluons. This should not be a serious drawback of the approach, 
having the r\ essentially an octet character under SU (3) transformations. 
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FIG. 1: The DA for the q (upper panel) and s (lower panel) flavor in the rj meson, at the initial scale Qq = 1 GeV 2 (dotted 
line) and after evolution to the scale Q 2 = 50 GeV 2 (full line). The asymptotic behavior is also shown for comparison (dashed 
line). 



We fix now the values of Qq, C3 and M. The Qo scale is closely related to the choice of the Aqcd value. We fix a 
scale of Qo = 1 GeV, together with Aqcd — 0.226 GeV, in analogy with the previous analysis 0. A natural condition 
to be satisfied is continuity between the low virtuality description of the rjTFF and the high virtuality description, 
provided by Eq (fT9]l . To minimize the model dependence, we use the parameterization of the CLEO collaboration 
|36j for the description of the TFF in the LV region: 



F 



LV 



77 v 



(Q 2 ) = 



647rr(?7 — > 77) 
(47ra) 2 TO 3 



1 



1 - -X- 



F(0) 

+ A?, 



(29) 



with F(0) = 0.272 ± 0.007 GeV" 1 , obtained using T(r) -> 77) = 0.510 ± 0.026 10" 3 MeV as given by the Particle 
Data Group 0] together with m„ = 547.85 MeV, and = 774 ± 29 MeV [Hj]. 

The value of the mass M can be obtained equating the rjTFF given by Eq. (fl"9|) at Q 2 — Qq, using as $,,(x, Qq) 
the one provided by the NJL model, to the value given, at the same scale, by the monopole parametrization Eq. ()I9|) . 



(Qo) 



Q 2 qF(0) 

+ A? 



dx 



x + 



M 2 



C 3 



(30) 



Finally, the only unknown is C3, for which several reasonable values have been used, as discussed in the following 
section. 



IV. RESULTS AND DISCUSSION 



We present now our results of the calculation of the rjTFF. The starting point is the r/DA, evaluated at the low 
energy scale of the model. According to Eq. (j2~4"]l . all we need is the value of the rj mass, the quark masses and 



7 




FIG. 2: Calculation of the transition form factor via the rjDA with M = 0.557 GeV, C 3 = 2.04 10" 2 GeV 3 and using 
Qq = 1 GeV 2 (full-line) compared with the available experimental data 0, H(| . The gray region describes the indeterminacy 
on C3. The dashed line represents the result obtained taking C3 = 0. The dotted line corresponds to the asymptotic value for 
the form factor. 




FIG. 3: The same as in Fig [2j but in the low virtuality region 
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Set I 


171.2 


430 


-186 


497 


541 


1157 


1.07 


0.603 


-0.664 


Set II 


184.2 


435 


-184 


515 


554 


1148 


1.07 


0.832 


-0.840 


Exp. 






-194 


495 


548 


958 


1.18 


0.828 


-0.848 



TABLE I: We show results for two different parametrizations of the NJL model for several physical quantities, together with 
their experimental or phenomenological values. Explicit expressions for these are reported in the Appendix. The masses and 
the quark condensate are given in MeV. 

the regularization parameter A. We use for the r\ mass the experimental value, m, n — 548 MeV . The quark masses 
m 9 ( s ) and the regularization parameter A have to be fixed within the NJL model. It is important to work in the 
Pauli-Villars regularization scheme, in order to preserve gauge invariance. Unfortunately, to our knowledge, all the 
available fits for the NJL model in SU(3) are done within the cutoff regularization scheme. The only exception is 
the paper by Bernard and Vautherin [45j . where anyway an approximate expression for the 7a (TOi, m.j,g 2 ) integral is 
used. Therefore, we have performed a new fit of the model parameters. The SU(3) NJL model gives a very good 
description for the meson properties |49| . but one has to be careful, since it does not include confinement. To avoid 
problems, we impose a value of m u > m v /2. The details of the model (whose Lagrangian is given by Eq. (jAip ) are 
given in the Appendix. Here it is worth to recall only that the model has five parameters, which can be chosen as the 
current quark masses, fj, u and /z s , the dressed quark masses, m u and m s , and the cutoff parameter, A. Our strategy 
for the fits has been: i) the m u mass has been fixed to 275 MeV; ii) fi u and A are determined by fitting and 
fir', Hi) /i s and m s are chosen looking for an overall good description of the strange sector. In table U two different 
sets of the relevant quantities, obtained by the above described fitting procedure, are reported. It is seen that their 
experimental values are reproduced very well. In the Set I, we have imposed the additional condition for the resulting 
eta mass: < 2m u . In this set, a very good description of masses in the strange sector is obtained, but paying 
the price of a worse description of f,f' s . In Set II, the description of the masses is slightly worse, but the f^' s are 
very well reproduced. One should notice that, in going from SU(2) to SU (3) within the NJL model, the number 
of parameters moves from 3 to 5. One has therefore at hand two more parameters for explaining three new masses, 
three new decay constants and a new quark condensate. Actually, for calculating the rjDA only m u and A, which 
have the same value in both sets, are needed, together with m s , which changes from 430 MeV in Set I to 435 MeV in 
Set II. With respect to these quantities, the predictions of the model arc therefore reasonably stable. The results for 
the rjDA are presented for the following values of the parameters: A = 740 MeV, m u — 275 MeV and m s — 435 MeV. 
It may be useful to reiterate that the obtained values of f^' s are not relevant for the present calculation, because we 
used the experimental ones. 

In Fig. [TJ the rjDA are shown. We observe that <^ (x, Qq) is peaked around the central point x = 0.5 while (/>* (x) 
is relatively flat. This is a consequence of the masses of quarks u and d, which are close to half the mass of the eta, 
while it is not the case for the mass of the strange quark. What is clearly seen is the following reasonable feature: the 
less bound is a system, the more narrow is its DA around the point x = 0.5. 

Our DAs have an infinite expansion in terms of the Gegenbauer polynomials. The firsts coefficients a^ s , defined 
in Eq. ([28]). at Q% = 1 GeV 2 are: 

a\ = 0.134 a s 2 = 0.377 

a\ = 0.352 a% = 0.245 (31) 

The a* coefficients are close to the values predicted by a fiat distribution. On the other hand side, we observe that 
a\ < a\, at variance with what is commonly used in the field. This feature is due to the narrow structure of 4>^ (x) . 
We can compare our results with the values used by other authors. The a\ and a| coefficients are to be compared to 
the parameter B — 0.3 used in Rcf. [24]. 

In Ref. [Tl| the values a\ = -0.06 ± 0.06, a\ = -0.07 ± 0.04 and a£= 0.22 ± 0.06 are given, but at Q 2 = 4 GeV 2 . 
From Eq. and using ^ = (5 §\ + <£*)/6 and ^ = (5 - 2 0«)/3 we obtain 

a\ (1 GeV 2 ) = 0.17 a% (l GeV 2 ) = -0.028 . (32) 

Evolving these results we have 

4 (4 GeV 2 ) = 0.14 al (4 GeV 2 ) = -0.022 . (33) 

The value for a| is therefore consistent with that used in Ref. [11], while some difference is found for the value of 
a\. In the pion case, if we use a flat distribution at Q\ — 1 GeV 2 , a value (4 GeV 2 ) = 0.31 is obtained. In Ref. 



9 



[25j it has been noted that the values for these parameters found in [ll[ suggest a very large SU(3) breaking between 
the DA of the ir° and the one of rj s , and a very little U(l) symmetry breaking between 77 s and r] 1 . Our results show 
the same structure of those of Ref [1JJ, at least for and for the big difference between a J and of. The origin of 
this difference is in the small value of a\ due to the narrow structure of <f>i(x), originated by the fact that m u is close 
to m v /2. At the same time, a small value of a\ explains a small value for a\ . One should remember that the present 
scheme reproduces the SU(3)_f and U(3)f symmetry breaking in the pseudoscalar meson sector. 

The results of Eq. (|32p are also in good agreement with those from [23|,[5(|. In the latter references, the authors 
give their results for Qq — 1 GeV 2 in terms of the quantities B\ and B\ , defined in [50j and related to our expressions 
as follows : a\ (l GeV 2 ) = -B\ + Bf/102 and a% (l GeV 2 ) = -B\, Using the numerical results of Ref [Hj], we 
have a\ (l GeV 2 ) = 0.149 ± 0.048 and a% (l GeV 2 ) = -0.0425 ± 0.0175, to be compared to our results, Eq. ([32]). In 
these papers, the coupling of a two gluon state to the singlet qq component of the rj mesons is introduced explicitely, 
providing a contribution, B%, which is an important part of the final result. In absence of gluons, the symmetry U(l) 
is not broken. In our case, the U(l) symmetry is broken through the 't Hooft interaction term [5l| introduced in the 
Lagrangian (|A1|) . making our results consistent with those of refs. (23l. [Hoj. 

As stated at the end of the previous section, once the r/DA has been obtained at the scale Q\ and evolved to Q 2 
according to Eqs. (|26 p -(|28 p . the only remaining unknown for the evaluation of the rjTFF according to Eq. (fH)]) is 
the constant C 3 of the higher twist term. To this aim, three different scenarios have been considered, corresponding 
to a contribution from this term to the form factor at Ql = 1 GeV 2 of 10% (C 3 = 1.02 • 10~ 2 GeV 3 ), 20% (C 3 = 
2.04 • 10~ 2 GcV 3 ) and 30% (C* 3 = 3.06 • 10~ 2 GeV 3 ). The cutoff parameter M varies between 487 MeV, for C 3 = 
1.02 • 10~ 2 GeV 3 , 557 Mev, for C* 3 = 2.04 • 10~ 2 GeV 3 , and 638 MeV, C 3 = 3.06 • 10~ 2 GeV 3 . We show in Fig H the 
obtained result for rjTFF and in Fig. [3] a detail of the region between Q 2 = and Q 2 = 10 GeV 2 . 

The results for the rjTFF, shown in Figs. [5] and [3] exhibit a very good description of the experimental data in 
the whole kinematic region. For completeness we have included in the figures the C 3 = case and the asymptotic 
vale for the rjTFF, Q 2 F^ Tv (Q 2 ) = (5^/2/^ + 2/ J( s ) = 0.181 GeV. It may be interesting to notice that the value 

of the asymptotic rjTFF is very close to the value of the asymptotic nTFF, Q 2 F 1 ~ I * 7T (Q 2 ) — V^fw = 0.185 GeV. 
It is also clear that, at variance with what happens in the tt case, to explain the eta data some C 3 contribution 
may be needed only in the region around Q 2 = 20 — 30 GeV 2 . Anyway, a complete discussion is obtained only by 
comparing the present results for the ttTFF with those of Ref. |8j. In the tt case, the C 3 contribution was crucial 
to reproduce the data in the region Q 2 = 10 — 20 GeV 2 and the calculated ttTFF crossed the asymptotic curve quite 
early (around Q 2 = 10 GeV 2 ) and with a significative slope. In the rj case, the situation is less dramatic: the higher 
twist term improves the TFF description only slowly, and the theoretical result crosses softly the asymptotic value 
around Q 2 = 40 GeV 2 . 

Another interesting point is the stability of the parameters. In calculating the rjTFF, we have adopted a procedure 
independent from that used in Ref. [8] , namely, C 3 and M have been fitted using the rj data only. The parameters used 
in both calculations have been Aqcd = 0.226 GeV and Qq — 1 GeV . Otherwise, in Ref Q, a fully model independent 
calculation was performed, choosing 4> v (x) = 1 on the basis of chiral symmetry. Here one is forced to choose a model 
for the description of the r\DA at Qq, and C 3 and M have been fixed within this model, independently from the tt 
case. Despite of this, the result obtained in the two calculations are quite consistent. Varying the weight of the higher 
twist term from 10% to 30% produces a change in C 3 from 1.02 • 10~ 2 GeV 3 to 3.06 • 10~ 2 GeV 3 in the rj case, to be 
compared to a variation of C 3 from 0.99 • 10 -2 GeV 3 to 2.98 • 10 -2 GeV 3 in the tt case. The agreement is impressive. 

On the other hand side, we found for the mass parameter a wider variation. In the rj case one gets M — 560 ± 
70 MeV, taking into account the uncertainty in C 3 , to be compared to M = 620 ± 70 MeV for the tt case. Despite of 
these differences, the results can be considered perfectly consistent with each other. The difference in the central value 
of M could imply that, for the pion, a larger contribution from the transverse momentum is expected with respect to 
that for the eta particle. It is indeed what has been obtained in Refs. [n], |24j- The values of M could be compared 
with the value of (k±) given by P. Kroll [11], (k±) ~ 710 MeV for the pion and (k±) ~ 390 - 440 MeV for the eta. It 
can be also compared with the /3 T parameter used by [24[ , which is related with the width of the gaussian distribution 
of transverse momentum used by these authors, with the values (3^ = 668 MeV for the it-quark and /?„• ~ 530 MeV 
for the s-quark. A comparison of our parameters, based on a quark-flavor decomposition of the relevant quantities 
(DAs, decay constants), with those used in Ref. [221 ] . obtained within a sing let-octet decomposition, is instead rather 
involved. The spirit of the present calculation and those of Refs. [22|, |24| are rather different. In our calculation, 
the known QCD evolution of the DA governs the Q 2 dependence of the form factor. The same Q 2 dependence is 
obtained, in Refs. [12, [Hj], through the kr dependence assumed for the light-cone wave function of the mesons. It 
is therefore significant that the two approaches provide similar results, describing probably, using different tools, a 
similar mechanism. 

In the present discussion, the result on 7* (q) — > P7, reported by BaBar for a time like q 2 — 112 GeV 2 , 
q 2 F 11 , 11 (q 2 ) = 0.229 ± 0.030 ± 0.008 GeV, has not been included. The reason is that the kinematics and dynam- 
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ics of this process is different from the ones studied here. There is no symmetry relating F^ ll -p(q 2 ) at one point q 2 
to F 77 »p(<5 2 ) in the point Q 2 = —q 2 - The coincidence is in the asymptotic value, which has been predicted for this 
process to be[Z(| —q 2 F 11 *p (g 2 ) = \/2/p (1 — 5a s (g 2 ) /3n), when the contribution coming from the a s (g 2 ) term 
could be disregarded. In the present scheme we obtain Q 2 F 11 *p(Q 2 ) = 0.19 GeV at Q 2 = 112 GeV 2 , which implies 
a very slowly growing behavior of the TFF even for these high values of the virtuality. 

In closing this section, it is useful to list items that prevent from using the same formalism for the description of 
the rfTFF. First of all, as it has been previously noted, the NJL does not include confinement. Therefore, if one uses 
the same expression, Eq. (1241) . in the rf case, an imaginary part will appear in the DA at some value of x. Secondly, 
the rf is basically a singlet state and it can mix strongly with the two gluons state or, later, with some cc component. 
These ingredients are not included in the present formalism. 

V. CONCLUSIONS 

In this paper, the rjTFF has been discussed in a formalism which connects a low energy description of the hadron 
involved with a high energy description based on a QCD perturbative formulation. The two descriptions are matched 
at some scale Q 2 . The scheme has been applied to describe the parton and generalized parton distributions with 
notable success |28l - l30l . l32l [35| and, in particular, to the nTFF in Q. The formalism selects therefore two regions 
of virtuality, separated at Qq. For Q 2 < Qq, use has been made of the experimental parametrization of the rjTFF 
data. This has been done to avoid model dependence in this region. At Q 2 > Qq, use has been made of a high 
virtuality description, which incorporates the following important physical ingredients: i) a rjDA obtained in the NJL 
model ; ii) a mass cut-off in the definition of the rfTFF from the rjDA, M, |3{ which, interpreted from the point of 
view of constituent models, takes into account the constituent mass, transverse momentum effects and also higher 
twist effects; in) an additional higher twist term into the definition of the rjTFF in the high virtuality description, 
parameterized by a unique constant, C3; iv) the two descriptions have to match at a virtuality Qq, a scale which is 
universal and should be the same for all observables. 

In section II it has been shown that the dominant, twist two, expression for the pseudoscalar TFF, given in Eq. 
(1141) has to be corrected, for including higher twist effects. The minimal correction would be the one given in Eq. 
(HI). 

The rjDA at Qq has been obtained in the NJL model. For that, the parameters have been adjusted for a good repro- 
duction of the rj sector with the Pauli-Villars regularization. The obtained fits represent an overall good description of 
the strange sector, not only for the masses, but also for the meson decay constants. It is worth to strees that, in going 
from SU(2) to SU(3), the number of parameters is increased by two, while the number of new physical quantities, 
included in Table [I] are seven. The obtained DAs in this model show consystency with other analyses, where the DAs 
are parametrized [111 [23L |24T ] . 

Using Qo = 1 GeV as matching point, the higher virtuality results of the rjTFF are well reproduced. The C3 term 
turns out to be relatively small. Its effect is to reduce the value of the contribution to the twist two rjTFF only for 
Q 2 < 5 GeV 2 . Its value, C3 = 2.04- 10 -2 GeV 3 for a 20% of higher twist contamination at Q 2 , is in perfect agreement 
with the one obtained for the irTFF in Ref. 8], C* 3 = 1.98 • 10~ 2 GeV 3 . Mor eover, the results are very stable with 
respect to variations of this parameter. 

The value obtained for M = 560 MeV is comparable with that of the pion case in Ref @, M = 620 MeV. The 
relative high value of M in both cases can be understood thinking that it includes the constituent quark mass, the 
mean value of the transverse quark momentum and other higher twist contributions. In turn, the higher value of M 
for the 7r than for the rj can be related to the fact that the contribution of the quark transverse momentum in the 7r 
case is expected to be more important with respect to the rj case [ill, H3] • 

The calculation proves that all the BaBar results can be accommodated in the present scheme, which only uses 
standard QCD ingredients and low virtuality data. It must be emphasized that, in order to have a good description 
for both 7r and rj, higher twist effects are important, as the modification from Eq. (|14l) to Eq. (jTSJ) signals. It 
must be also noted that the matching scale is as high as 1 GeV, a feature already found in the description of parton 
distributions when precision was to be attained. With these ingredients, the calculation shows an excellent agreement 
with the data. 

Let us conclude by stressing that we have justified the formalism developed in Ref. [§[ to describe the irTFF and 
we have extended it to the rjTFF. The idea of the approach is that one can use models or effective theories to describe 
the non perturbative sector, and QCD to describe the perturbative one. In here, we have preferred to use data for the 
low virtuality sector to avoid model dependence, but in building the rjDA at Qq we have used the NJL model. Higher 
twist effects (parametrized in our case by M and C3) are small but crucial in order to attain an excellent description 
of the 7r and rj experimental results. 
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Appendix A: The SU(3) NJL model for pseudoscalars mesons. 



In calculating the T]DA, the minimal extension of the NJL model for describing pseudoscalar mesons in SU(3), 
proposed in Ref [HI, has been used: 



8 

C = q (x) (iQ-n) q (x)+Gj2 [m a l) 2 + (Qi^X a qf 



a=0 



K [det(g(l + 75 )g)+det (g(l- 7s ) g )] 



(Al) 



where fi = diag[/i u ,/Xd,/z s ] is the matrix of the current quark masses and A a , a = 0, ...,8, are the SU(3) generators. 
SU(2) will be considered a good symmetry, and, therefore, fi u = As it is well known, the first consequence of the 
scalar interaction term is to provide the constituent quark masses, m u = m<j, m s , different from the current ones. 
The main results are summarized here, while the reader is referred to the section IV-B of Ref. [38| for details. 
By defining the integrals: 



I\ (m) = i 
I 2 {rrii.m^q 2 ) = i 



d A p 



1 



(27r) 4 (p 2 -m 2 + ie) 



d 4 p 



(2tt) (p2_ m 2 +2e) (p-qf- 



\2 2 



(A2) 



the constituent masses are given by 

m u = fi u + 16 G N c m u h (m u ) + 32 K N 2 m u h (m u ) m s h (m s ) , 

m s = Ms + 16 G N c m s h (to s ) + 32 K N 2 [m u h (m„)] 2 , (A3) 

where N c is the number of colors. The vacuum expectation values for the condesates of the quarks of flavor qi are 

((Mi) = 4TOjIi (mj) - (fit) , (A4) 

where the expectation value of the qq in the perturbative vacuum has been substracted from the expectation value 
in the true vacuum. The last term in Eq. (|A4[) is negligible in the u — d quark sector, while it becomes important in 
the strange sector. 

The next step is the description of the pseudoscalar states. For the pion and kaon case, by defining the quantities 



K 3 = G + 2N c Km s Ii(m s ) , 
K G = G + 2N c Km u h (m u ) , 



and 



D {mi,mj,q 2 ) = Ii(mi) + hirnj) + (^(rni - mj) 2 - q 2 ^j I 2 {ra il m,j,q 2 ') , 

one has that the pion and kaon masses are obtained solving the equations: 

l-8N c K 3 D(m u ,m u ,ml) =0 , 
1-8N C K 6 D (m„, m s , m 2 K ) = . 

The couplings of the pion and the kaon to the quarks are given by: 



(A5) 



(A6) 



(A7) 



2 

9nqq 



d , o 



2 

9Kqq 



4A^ C -j~p D (m u ,m s ,q 2 ) 



(A8) 



and the decay constants are: 



F x = —iN c g^ qq m u I 2 (m u ,m u ,ml) 



(A9) 
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Fk — 2Nc9 2 Kqq |( ms + m u ) D (m U7 m s , m 2 K ) -2m u h {m u ) - 2m s I 1 (m s )} , (A10) 
m K 

where F„,k = U,k/V%- 

The 77 particle deserves a more careful discussion, due to its mixing with the rf particle. Working in the flavor basis, 
one can define 

K q = G-2N c Km s h(m s ) , 

K qs = -2V2 N c Km u h(m u ) , (All) 

JC = K q G — K qg . 

The interaction in the r\ — rj sector can be described by the expression 

( l75 A*) [My] (ij 5 X j ) (A12) 
with i,j = q, s, X q = diag[l, 1, 0], A s = diag[0, 0, y/2] and the interaction matrix is given by 



1 



[Ma] = -jj- I aqq °" I , (A13) 
U v V a qs a ss I 



with 



a qq = 2(K q -81CN c D(m s ,m s ,q 2 )) , 
a ss = 2(G-8!CN c D(m u ,m u ,q 2 )) , 



qs -2K qs , (A14) 



D v (l 2 ) = ( a ii a ss - a 2 qs )/(4IC) . 
The r\ mass is obtained solving the equation 

D v (m*) = . (A15) 
In a neighborhood of q 2 = m 2 the interaction can be written as 

(- sin X s + cos 4> X q ) ~ 9mq (- sin <j> X s + cos <f> X q ) 



9 2 
q- — mf, 

-f (e v X s + X q ) (e ri X s + X q ) , (A16) 



a. 
'n 



with e v = a qs /a qq . In obtaining the right hand side of this equation, use has been made of Eq. (|A15j) . which implies 
a ss = a qs /a qq . From (|A16[) one has 



ff. 



, (1 + 4) a ii 



mi 



dD ri /dq 2 



(A17) 



For the flavor decay constants, one has 

F q = -12g mq costf)m u I 2 (m u ,m u ,m 2 ) , 

= 12g vqq sm(f) m s h (TO s ,m s ,m^) , (A18) 

where F«> s = f q ' s /V2. 

We need to evaluate the integrals defined in Eq. (IA2[) . Due to the point-like character of the interaction, the 
lagrangian Eq. (|A1[) is not renormalizable and a regularization procedure for these integrals has to be defined. We 
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use the Pauli-Villars regularization in order to render the occurring integrals finite. This means that, for integrals 
like the ones defined in Eq. (|A2|) . we make the following replacement, 



h (tth) — > X! Ci Jl 12 (mi,mj,q 2 ) — > ^ c e I 2 (Me ti , M^j, q 2 



(A19) 



£=0 



with M%j = m 2 + £A 2 , Co = c 2 = 1, c\ = —2. Here, for simplicity, we choose the same A value for the strange and 
the nonstrange sector. According to these prescriptions one finds 



i 



16 7T 2 



M 2 „ M 2 

-2 M 2 i In — ii + M 2 , In — 3i 

mf ' mi 



(A20) 



1 2 

him^m^q 2 ) = — ^ c e 

i=0 



In 



In 



III; 



M 2 



hi 



M, 



1,3 



M? 

lvl l,i 



(A21) 



with 



-X[Mf ti ,Mfj,r 



arctan ■ 



Mf. 



Ml, 



+ arctan ■ 



Mh + Ml, 



^j-\(Ml,Ml p q>) ^j-X^Ml^) 



(A22) 



where A 

(Mf,i, M l^Q 2 ) is the Kallen lambda. 
Now, we fix the parameters of the model. Looking at the lagrangian, we have a five parameters model, fi u , fi s , G, 
K and A. Nevertheless, it is more intuitive to organize the fit of the parameters in terms of fj, s , m u , m s and A, 
using equations (|A3j) to determine G and K. We impose m u — 275 MeV, in order to have m® xp < 2m u . Then, [i u and 
A are obtained in recovering the values of and m^. At this step one has 

\3 



6.69 MeV, m u = 275 MeV, A = 740 MeV 



138 MeV, F„=92.2 MeV, 



(-227 MeV)- 



determining the SU(2) sector. Then, /z s and m s have been fixed by requiring a good overall fit of masses (ma, 
m v , m,f ) and decay constants (Fk, F^, F"). In table [J two different sets of parameters are given, together with the 
obtained results. Using set I, by imposing m v < 2m u , a good agreement for the masses and a slightly less good 
agreement for the F^ ,s is obtained. On the other hand, in set II a very good agreement for F^' s is obtained, with a 
slightly worse result for the masses. 

In the light-front calculation one needs the integral 



h i,m„m p 



5 (x - 1 + 



d 4 fc 

(2tt) 4 [(k- P) 2 - m 2 + ie] (k 2 - mf + it) 



(x)0{l - x) 



(Air) 2 ^ 
y ' l=i 



ce In 



mf — (1 — x) x m 



M% i - (1 -x)xm 2 P 



(A23) 



Clearly, all one needs to calculate the rjDA are the 77 and quark masses and the value of the cutoff parameter. For 
the DA calculations, the values m u = 275 MeV, m s — 435 MeV, A = 740 MeV and m n — 548 MeV have been chosen. 
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